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Abstract Single and double torsion oscillators have been used to measure the 
anomalous change in resonant frequency and accompanying dissipation in solid 
O ', ^He. We present a glass description of the mechanical anomalies found in tor- 

^ • sion oscillator measurements. Our results show that it is not necessary to invoke 

a supersolid interpretation to explain these mechanical anomalies. Previously, we 
demonstrated that the back-action of a glassy subsystem present in solid ^He can 
account for frequency change and dissipation peak in single torsion oscillator ex- 
periments. Here, we show that the same glassy back-action can explain the experi- 
C^ • mental results of the composite torsion oscillator developed by the Rutgers group, 

which measures the response of solid *He at the in-phase mode /i = 496 Hz and 
out-of -phase mode /2 = 1 173 Hz. 

C3 ■ PACS 67.80.B-, 67.80.bd, 64.70.P- 

fi I Keywords solid helium, supersolid, glass, torsion oscillator 



1 Introduction 



(N 
> 

Qs I The controversy over the origin of the anomalous low-temperature signature in 

Q>^ . solid ^He continues since it was first discovered in ultrasound measurements by 

Goodkind and collaborators-^^ and eventually detected in torsion oscillator (TO) 
experiments by Kim and Chan.'*-^ In earlier work, we argued that the origin of 
f~^ ■ the TO anomaly is not that of a non-classical rotational inertia effect (NCRI), but 

f^ I rather due to a mechanical effect. We examined the response of a nonuniform solid 
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with defects leading to transient elastic dynamics at high temperatures. i'i' 
These defects are quenched at lower temperatures leading to increased mechanical 
stiffness and hence increased rotational speed of the TO. Furthermore, a nonuni- 
form solid with a small concentration of defects can explain the maximal damping 
at the crossover temperature where the putative NCRI was reported. Our defect- 
based scenario is quite simple. At high temperatures, the transient dynamics is 
very rapid relative to the time scale set by the TO, while at temperatures far lower 
than the putative supersolid transition, the transient dynamics is much longer. As 
in classical damped oscillators, when the transient dynamics due to internal dissi- 
pation matches the impedance of the external driving force, the dissipation peaks. 

A "freezing" between a high-temperature dynamic state and a low-temperature 
immobile state may be characterized by either (i) a discontinuous first order or 
continuous transition in annealed systems with no disorder or (ii) an observed dy- 
namic crossover in amorphous systems with disorder (spin-glass) or self-generated 
quenched disorder due to rapid cooling (as in structural glasses). Freezing of rapid 
degrees of freedom can be triggered by a distribution of activation barriers; as we 
will briefly review below, activated dynamics corresponds to a particular choice of 
the relaxation time for glasses. Such a distribution of barriers also lies at the heart 
of theories describing insulating "electron glasses" in which electrons are pinned 
by disorder. ' ^ Some theories of structural glasses assume that there is no true glass 
"transition", but rather a dynamic crossover at which the relaxation time becomes 
long relative to experimental time scales. Yet other theories assume that an ideal 
glass exists in the thermodynamic limit but that, experimentally, it is largely inac- 
cessible by divergent equilibration times. Notwithstanding the lack of consensus 
on the origin of universal slow dynamics in numerous amorphous systems, the in- 
crease in observed relaxation times is due to pinning or freezing of some form of 
degrees of freedom that are dynamic at high temperature. ^^ 

Several theories of glasses assume that elastic defects can entangle and become 
pinned and consequently form a sluggish state at low temperatures.'^ An early 
theory of defect dynamics investigated the consequences of a low-temperature 
pinning of dislocations. '^ Although, as emphasized above, the underlying mecha- 
nism of low temperature dynamics in these amorphous systems is not universally 
agreed upon, their empirical behavior has precise universal features,— which we 
employ here. Salient features of amorphous systems that are pertinent to our study 
include: 

• A nearly universal Vogel-Fulcher-Tammann (VFT) form describing the in- 
crease in relaxation time of glasses as temperature is lowered. The VFT form has 
been used to fit data collected from many glasses, although it has no unique the- 
oretical basis. Specifically, the VFT expression T = Toexp(zi/(r — 7b)) includes 
activated (Arrhenius) dynamics (Tq = 0) dominated by energy barriers A. Thus, 
the general VFT analysis that we employ includes the specific case of such ac- 
tivated dynamics; this particular case was discussed in our earlier fits to the TO 
data.^ Indeed, activated dynamics triggered by pinning of defects (dislocations) 
may occur in solid ^He. '^ Similarly, increased dynamics relative to an activated 
form (To < 0) suggests an increase in motion due to quantum fluctuations. 

• The activated dynamics in amorphous systems is characterized by a broad 
distribution of relaxation times. The empirical response functions, which capture 



this distribution, are consequently of the Davidson-Cole or Cole-Cole form which 
we invoke in our analysis. 

• Omni-present memory effects and hysteresis are particularly important in 
these systems and characterize the observed dynamic crossover including "aging" 
behavior.— li^ii^ Such memory effects are of current technological importance.— 

• Finally, low-temperature specific heat contributions that scale linearly with 
temperature in insulators have been key indicators for glassiness.^''^^ 

It is important to emphasize that our approach assumes that only a small 
glass-like fraction present in solid "^He exhibits a freezing transition of a high- 
temperature mobile component with a distribution of relaxation times. ^ Since we 
first proposed a nonsupersolid interpretation, several other disorder-based theories 
have been put forward to explain the large dissipation seen in torsion oscillator or 
shear modulus experiments. Yoo and Dorsey proposed that solid '*He displays vis- 
coelastic properties. ^^ Syshchenko et al.^^i^ and Iwaza^^ argued that the shear 
modulus measurements can be understood in terms of dislocation motion and 
the formation of a low-temperature dislocation network based on the theory by 
Granato and Liicke. ^^ All these models have in common disorder, but do not re- 
quire glassiness. Similarly, Syshchenko ^^ speculated that the observed dissipation 
peak is associated with boiling off of ^He atoms from dislocation lines. The ^He 
unbinding causes a depinning with increasing temperature that generates a dis- 
tribution of mobile dislocation lines responsible for softening of the shear mod- 
ulus and increased dissipation. Very recently, Reppy^^ reported a nonsupersolid 
scenario in TO experiments, where induced disorder affected primarily the high- 
temperature response above the putative supersolid transition. A result consistent 
with the dynamics of highly disordered or glassy systems. So far attempts to of- 
fer a unified picture within elasticity theory of uniform solids, i.e., how applied 
shear stresses affect torsion oscillators and elastic coefficient measurements, have 
been unsuccessful. 24i29j30 -pj^^ j^^^y ^^ j^g ^^ ^ neglect of anelastic contributions 
in these approaches. 

Since a true phase transition is a bulk effect of the sample, we investigated 
previously thermodynamic signatures as reported by specific heat and pressure 
measurements 3 1 ,32,33,34,35 q^j. thermodynamic analysis in terms of two-level sys- 
tems is consistent with a glassy interpretation of the linear-T dependence in the 
specific heat below ^ 100 mK and the quadratic-T dependence in pressure mea- 
surements of the equation of state, P{T), of an otherwise perfect Debye solid, ^^i^^ 
A Debye solid exhibits a specific heat C '^T^ and a pressure dependence P oc T^. 
It is notable that the associated excess entropy is at least three orders of magnitude 
too small to explain an NCRI effect of a couple of percent in the TO experiments if 
it is caused by uniform Bose-Einstein condensation. ^^-^^ Finally, we demonstrated 
previously that single torsion oscillators can be described by using a distribution 
of glassy relaxation times, for example, the well-known Cole-Cole distribution. 
In Fig. [U we show excellent agreement between model calculations for a glassy 
back-action and the data by Hunt and coworkers for a single TO.^^ Details of the 
model have been described in Ref. ^ The important message of this comparison is 
that 100% of the TO signal can be ascribed to a glassy back-action response of 
mechanical origin with no need for a supersolid component. 

Here we focus on the TO experiments reported by Aoki and coworkers.—"^ 
The unique design of two coupled torsion oscillators has opened up a window in 



574.48 




Hunt et al. 

Cole-Cole distribution (a=1.85) 



■e- 



-e G- 



100 200 

T(mK) 



10 



O 



300 



Fig. 1 (Color online) Torsion oscillator resonant frequency (black, left axis) and dissipation 
(blue, right axis) vs. temperature. The experimental data (symbols) for a single torsion oscillator 
are from Hunt et al.^^ Details of the glass model (curves), the Cole-Cole distribution function 
of Debye relaxors with exponent a = 1.85, activation banier A = 408 mK, and other model 
parameters are given in Ref — 



frequency space to explore the dynamic response of solid '^He at two very different 
frequencies. However, these measurements have been difficult to explain within a 
supersolid or glassy scenario alike. In the past, these difficulties have been used 
as arguments against a glassy interpretation. 37,38^29 j^^ composite design allowed 
for the measurement of an in-phase, /i = 496 Hz, and out-of -phase, /2 = 1173 
Hz, resonant frequency on the same solid ^He sample. Surprisingly, the observed 
changes in frequency and dissipation as a function of temperature showed that the 
relative magnitude in frequency shift, Afi/fi = (/;(0K) -/;(0.3K))//,(0.3K), is 
nearly frequency independent. At the same time, the dissipation peak shifted only 
slightly to higher temperatures at the higher frequency. These observations are 
contrary to the standard theory of supersolidity, with the hallmark of dissipation- 
less superflow of vacancies or interstitials. Even when invoking a vortex picture 
for supersolidity to account for the anomalously large dissipation, it fails to ac- 
count for the unchanged relative shift in resonant frequency Afi/fi w 4/2//2, 
which is a direct measure of the NCRI. Any theory of superfluidity predicts that 
the superfluid fraction is more rapidly destroyed when perturbed at higher fre- 
quency. This is certainly at odds with the experiments by Aoki and coworkers and 
deserves more attention. Simply extrapolating results of the single TO to higher 
frequencies fails to explain the experiments. The puzzle is how to explain the out- 
come of the composite double TO experiment. 

In this paper, we present an interpretation of these measurements using a 
glassy back-action term in the response function of the double TO similar to our 




Fig. 2 (Color online) Sketch of the double torsion oscillator modeled in Eq. (T). The upper 
moment of inertia (/| ) is the dummy bob, while the lower moment of inertia (h ) is the cylindrical 
pressure cell that can be loaded with '*He. The stiffness of the BeCu torsion rods is given by ki 
and k2 with ki fs ki by design. 



earlier work for the single TO. We provide a consistent interpretation and model- 
ing of the composite double TO within the glassy framework, where we postulate 
the presence of a glassy subsystem in solid ^He. We descibe the experiments by 
Aoki in terms of two coupled mechanical torsion oscillators, shown schematically 
in Fig. |2] Our model can explain the measurements, though requires anomalous 
damping and a strongly frequency dependent glass term. Note that anomalous 
damping is already required for the empty cell (no "^He present!) to explain that 
the quality factor Q2 of the high mode is lower than Qi of the low mode at 0.3 
K, its origin may be due to anelasticity or non-rigid connections between the tor- 
sion rod or to a nonlinear feedback of the capacitive drive onto the dummy bob. 
Additionally, our results point toward '*He slipping at the container walls of the 
cylindrical pressure cell consistent with the very small, frequency-independent 
relative shifts Afi/fi w Af2/f2 « 1.7 • lO"*'. 



2 A model for the coupled double oscillator 

We model the coupled double torsion oscillator of the Rutgers group, sketched in 
Fig. 121 by the following system of equations for the torsion angles: 

/i0i(O + yi0i(O + fci0i (0+^2(01(0 -02(0) = ^(0, 

l202{t) + 7202(1) +k2(02{t)-0l(t)) = J dt' git- t')02it'), (1) 

where 0,(f ) are torsion angles, Yi are damping coefficients, k, are torsion rod stiff- 
nesses, g{t) is the glass back-action term, and F{t) is the applied external torque. 
The subindex "1" refers to the upper or dummy bob in the experiment, while "2" 
refers to the lower oscillator with the pressure cell that can be loaded with solid 
^He. After Fourier transformation of Eq. ([Tl, we obtain 

{~h(O^-iYl(O + ki+k2)0l{(O)-k202i(O) =F{(0), 



{-I2(o^-ir2(o + k2-g{(o))02{a))-k20i{(o) =0. (2) 

For a strongly underdamped oscillator and a small glassy back-action, it suffices 
to solve first for the bare resonant frequencies and later include perturbatively 
damping and glass terms. The bare resonant frequencies {F = 0) are 
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Next, we expand the system of equations in (|2]l around the bare resonant frequen- 
cies OJ by inserting ft), = (O^ + 4 ft), and solving the secular equations for 

_ gi(ft)f)/?i(ft)°) + /[g2K'')/?i(ft),") + ft),°(yi/?2«) + y2/?i(ft),°))] 

2ft)0[/i/?2(ft)f)+/2/?l«)] ' ^^ 

where we introduced ancillary functions for compactness of notation 

g{(0)=gi{(0)+ig2{(i)), (5) 

/?i(ft)) = -ft)2/i+fci+fe2, (6) 

R2{(0) = -aP-l2 + k2. (7) 

It is now straightforward to calculate the resonant frequencies and dissipation of 
the loaded oscillator 

f. - ^, 

At high temperatures, where ft)T — > and g2 ~ the dissipation is given by 

1^ yi/?2«) + y2^i(t»,°) 

^''-~ft)«[/l/?2«)+/2/?l(ft)°)]' ^ ^ 

Since the empty pressure cell shows significant temperature dependence in its 
response, we follow the experimental analysis and subtract the T-dependent empty 
cell response and correct for the differences due to the small difference in I2. 
Hence, we compare against corrected shifts relative to a high-temperature value, 
e.g., at 300 mK, 

Afi = fi{T)-fi{03K). (11) 

Within our notation we find that approximately QrJ_ w gr' (0.3 K). 

Our model of the coupled composite oscillators can describe the temperature 
behavior of the Rutgers experiment if we make further assumptions about the 
glassy back-action g(ft)). We model the dynamics of the glassy subsystem by 
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Fig. 3 (Color online) Frequency and dissipation in double torsion oscillator by Aoki et al. — 
(symbols) compared with glass theory (lines). Panel (a): Temperature dependence of resonant 
frequency shifts Afi (black, left axis) and Afj (blue, right axis). Panel (b): Temperature depen- 
dence of dissipation QJ^ (black, left axis) and Qj' (blue, right axis). The experimental data 
were corrected for the strongly temperature-dependent background of the empty cell— 



Here the glass term is ^(ft)) = go ( "ti ) and the Vogel-Fulcher-Tammann relax- 
ation time is t{T) — Toexp[zi/(r — Tq)] for T > Tq. Tq is the ideal glass transition 
temperature, which is below the temperature where the peak in dissipation occurs. 
The parameter A is an average potential barrier of the glass. Finally, at tempera- 
tures T <Tq the glassy subsystem freezes out and T becomes infinite. 

The parameters /, and fc, can be determined from the bare resonant frequen- 
cies ff = (O^ /2n in Eq. OJ. In addition, the damping coefficients Ji can be ex- 
tracted from the high-temperature dissipation 2rj_. Finally, the glassy back-action 
^(ft)) with its parameterization accounts through t{T) for the temperature depen- 
dence oi Afi and Q,~ ■ Our phenomenological glass theory, applied to the coupled 
double oscillator, explains both frequency shift and dissipation peak for in-phase 
and out-of -phase torsional response. ^^ Data for in-phase frequency f\ = 496 Hz 
and out-of -phase /2 = 1173 Hz are shown in Fig. [3] plotted against the tempera- 
ture. The data are found to be well described by our model of coupled oscillators 
presented in Eq. (|2| with model parameters: /i = 620.3 nNs^, I2 = 233.2 nNs^, 
h = 9.125 N, k2 = 8.355 N, /i = -16.93 nNs, 72 = 10.43 nNs, go = 15.60 
IXN, To = 4.906 ^s, A = 386.2 mK, To = -32.73 mK, a = 1.74, p = 1.77. 
The obtained values for moment of inertia and rod stiffness agree well with es- 



timates for moments of inertia and BeCu rods.— At the highest measured tem- 
peratures, where the glassy contribution is small, these parameters result in reso- 
nant frequencies /i (0.291 K) = 495.82813 Hz, /2(0.300K) = 1172.8159 Hz and 
dissipation Q^' (0.291 K) = 0.820-10-^, g^^ (0.300 K) = 2.051 • 10"^ in excel- 
lent agreement with experiment. The experimental values are /[^'"' (0.29 IK) = 
495.82811 Hz, /f'"- (0.300 K) = 1172.8158 Hz and (ej^^)^^'" (0.291 K) =0.806 • 
10-^ (e^^)'''^/"(0.300K) =2.061 ■ 10"^ 

It is worth pointing out that an anomalous damping coefficient y\ ^ — 72 is 
needed to explain the anomalous behavior of increased dissipation with increased 
frequency. Note that this anomalous damping is already required to describe the 
unloaded pressure cell, so it is unrelated to the properties of solid ^He. After load- 
ing the cell with solid ^He the dissipation ratio becomes Q^^ /Q~[^ = 2.5 at 300 
mK with frequency ratio /2//1 = 2.37. A negative To is indicative of strong quan- 
tum fluctuations leading to an avoided glass transition. This behavior resembles 
that of a paramagnetic system with a negative Curie -Weiss temperature in the pres- 
ence of antiferromagnetic fluctuations. It will be interesting to see if the zero-point 
motion of "^He is indeed responsible for this behavior. 

Finally, the comparison in Fig. [3] shows that an explicit frequency-dependent 

back-action must be used with ^(ft)) = go [-%] and p = \ .11 to account for the 

experimental fact of Af\/ f\ « Af2/ fi, i-c, the relative frequency shift remains 
unchainged with changing resonant frequency. In the past, when we studied sin- 
gle oscillators, it sufficed to parameterize ^{(o) by the single resonant frequency, 
namely, ^{(o) w ^((0°) = go. 

It will require further studies to sort out whether the negative damping is re- 
lated to an anelastic or non-rigid torsion rod or nonlinear feedback from the ex- 
ternal drive. Notice that the combined oscillator system is underdamped and dis- 
sipates energy with Q, > 0. On the other side, a frequency exponent of /? = 1 .77 
in the glass term ^(ft)) points toward slip of solid '^He at the cylindrical con- 
tainer wall. Theories describing solid ^He in torsion oscillators as a viscoelastic 
material ^^ or two-level systems moving through a solid matrix^^J^^ili predict an 
explicit frequency exponent of p =A for the back-action term, when assuming a 
no-slip boundary condition at the cylindrical container wall. So it is reasonable 
to expect that the slip of ^He at the container wall will lead to a reduction in the 
power of the frequency dependence of ^(w). 



3 Conclusions 

For the first time, we present a quantitative explanation of the anomalous fre- 
quency and dissipation dependence reported for the double torsion oscillator with 
cylindrical sample chamber by the Rutgers group invoking glass dynamics. So 
far these experiments have been at odds with supersolid and glassy interpreta- 
tions. The conventional single oscillator models fail to account simultaneously 
for the low and high resonant frequency data. Our studies of the coupled oscil- 
lator show that the observed shifts in resonant frequencies and dissipation are 
in agreement with a glassy back-action contribution in solid ^He if one includes 
anomalous damping in the dummy bob and explicit frequency dependence of the 



glassy back-action term. Surprisingly, already the double torsion oscillator with 
the empty pressure cell requires a negative damping coefficient for the dummy 
bob to accurately describe frequencies and dissipation at 300 mK. Clearly, more 
dynamic studies of solid '*He with coupled oscillators are needed to determine the 
frequency dependence over a larger frequency range and for different cell designs 
to identify the origin of the anomalous damping and the frequency dependence of 
the back-action term. It remains to be seen if a quantitative supersolid interpreta- 
tion of the double oscillator experiment is possible. 
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